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ABSTRACT 

We study the wave optics feature of the gravitational microlensing by a binary system com- 
posed of parent star and a planet. In the binary system, near the caustic lines multiple images 
play the role of secondary sources for the observer, in analogy to the double slit Young's ex- 
periment. In the case of having coherent wave fronts from the source on the lens plane, images 
can produce diffraction pattern on the observer plane. For the binary lensing system we have 
two modes of close and wide images around the planet and lens star and these images can pro- 
duce two different types of fringes with the high and low frequencies on the observer plane. 
By taking into account the finite size of the source star, enhancements in the diffraction fringes 
get dimmer. For the observational prospects, we study this effect for the SKA project in the 
case of resonance and the high magnification exoplanet channels. This method can partially 
break degeneracies between the lens parameters. 



1 INTRODUCTION 



Gravitational lensing addresses to the bending of light rays by the 
gravitational effect of foreground mass. Depending on the distribu- 
tion of mass on the lens plane and the relative location of the lens 
and source from the observer, multiple images or distortion of the 
background source can be formed. In the case of star-star lensing 
inside the Milky Way, the separation of the images are less than few 
milliarcsecond hence they are unresolvable for the ground based 
telescopes. This type of gravitational lensing is so-called gravita- 
tional microlensing. 

In 1936 Einstein derived the gravitational lensing equation, 
but it took several years until 1979 that first gravitational lens- 
ing was observed. The source of this lensing was a quasar an d 
observation was done in the radio frequencies jWalsh et al 19791) . 
A few years later Paczyski proposed studying MACHO 's Q pop- 
ulation in the Galactic halo by the gravitational microlensing 
dPaczvski 19861) . His suggestion was observing stars in the Large 
and Small Magellanic Clouds and by counting the number of mi- 
crolensing events and measuring the transit time (Einstein crossing 
time), it would be possible to obtain the contribution of MACHOs 
in the Galactic halo. Beside the dark matter studies, the other in- 
teresting astrophysical application of microlensing was suggested 
by Moa and Paczynski (1991), using gravitational microlensing for 
the discovery of exoplanets. 

The microlensing effect by single or multiple lenses was 
mainly studied in the geometric optics. An important study on the 
wave optics feature of the gravitational lensing was done by Oha- 
nian (1983) who investigated the amplification of a radio point 
source when a galaxy acts as a gravitational lens. He showed that 
wave optics smoothes singular features of the light curve in the 
caustic crossing. In another work, Jaroszyski and Paczyski (1995) 
studied the caustic crossing of object like Q2237+0305 by a galaxy 
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composed by the individual stars. Studying the diffraction images 
of this system could put limit on the size of source. In Galactic 
scales the observation of the wave optics features also has astro- 
physical application as studyin g the limb darkening o f the small 
size sources as white dwarfs i Zabel & Peterson 20031) . Recently 
Heyl (2010,201 la,201 lb) argued about the possibility of detection 
of wave optics signals by the individual sub-stellar objects in the 
microlensing light curve. 

In this work our aim is to extend the application of the wave 
optics to the conventional observation of extra solar planets by the 
gravitational microlensing. Here we take a binary system composed 
of a lensing star and planet. Crossing the caustic lines of this system 
by the source star, produces high magnification in the light curve. 
For some of the cases due to the small separation of the images on 
the lens plane, the gravitational lensing system resembles to a mul- 
tiple slit optical system in the astronomical scales. Having coherent 
condition for the wave fronts on the lens plane, the result would be 
a diffraction pattern on the observer plane. We study the applica- 
tion of this method both for the resonance and high magnification 
channels of the exoplanet observations. The observation of the en- 
hancements and transition of the fringes will enable us to break 
degeneracy of the lens parameters. We also study the prospect of 
observation of the wave optics signal by future Square Kilometer 
Array (SKA) project. 

In section (O, we used the wave optics formalism and calcu- 
late the light curve of a point like and the extended source, lensed 
by a binary system. In section 10 we perform semi-analytic ap- 
proach of the wave optics in the microlensing near the fold caus- 
tics and study the conditions having the wave optics features. In 
section Q we discuss about the possibility of detection of wave 
optics signature of the extra-solar planets in the radio and microm- 
eter emissions of the source stars and prospect the observation of 
the light curves by SKA. We also discuss on the physical parame- 
ters that we can get from the light curve and degeneracy breaking 
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in the resonance and high magnification channels. The conclusion 
and summary is given in section 10. 



2 WAVE OPTICS IN GRAVITATIONAL LENSING 

In the geometrical optics, the locations of the images can be ob- 
tained from solving the lens equation 



y = x-a(x), 



(1) 



where x is the angular location of the the image, y is the angular 
position of the source and a(x) is the deflection, given by 



a(x) = - 

7T 



k(x.')- 



-d x 



(2) 



where k(x) = E(a;)/E cr , E" 1 = (A-KGD d D da )/{c 2 D s ) and D a 
, Dds , Dd are the source, lens-source and lens distances, respec- 
tively. T,(x) is the surface mass density of the lens. Another way 
to get the position of the images from the position of the source is 
using the Fermat principle. For the stationary points of the Fermat 
potential, the position of the images obtain as 



V x <Kx,y) =0, 



where Fermat potential is given by 

<X x ,y) = ^(y- x ) 2 -V'(x), 

and deflection angle is 



q(x) = V x V(x). 



(3) 



(4) 



(5) 



Here x and y are the image and source positions on the lens plane, 
normalized to the projected Einstein radius in each plane. Using the 
determinant of Jacobian transformation J(xi), we can obtain the 
magnification of a point like source mapping from the lens plane as 
follows: 



H(xi) 



\detJ(xi) 



\<j}ll(Xi)(f>22(Xi) - (j)\ 2 {Xi)\ ' 



(6) 



where Xi is location of ith image. For an extended source, we can 
calculate the magnification over source area. There are advance nu- 
merical methods such as using Green's theorem where two dimen- 
sional integration over the source t ransfer to one d imensional inte- 
gration on the boundary of images l lDominik 2007h . 

In the wave optics, from the Huygens principle, every point on 
the lens plane can be considered as a secondary source and the am- 
plitude of the electromagnetic wave on each point of the observer 
plane is composed by the superposition of the infinitesimal sources 
on the lens plane. From the Kirchhoff integral we can obtain the 
amplitude of the electromagnetic wave F^ v , knowing the boundary 
condition on the lens plane and when the d istance to the observe r 
screen is large compare to the wave length, feorn and Wolf 20021) . 
Multiplying the overlapped electro-magnetic filed to its complex 
conjugate results in the magnification of the light on the observer 
plane as follows l lSchneider et al. 19921) : 



4tv 2 



(7) 



where /0(x, y) corresponds to the elapsed phase for the elec- 
tromagnetic wave emitted from source, deflected from x and re- 
ceived by the observer. / = 2kR s and R s = 2GMt/c 2 is the 
Schwarzchild radius of the total masses. This formula is obtained 
for a monochromatic electromagnetic wave. 



In order to show the compatibility of the wave optics with the 
geometric optics and also transition from the wave optics to the 
geometric optics, we do Taylor expansion of the Fermat potential 
on the lens plane. For a given location of the source, the expansion 
around the image is given by 



(x,y) = 

+ 



,«>, 



+ (x - Xi).V a 



^ (0, + i[(-i 



{x 2 - X 2 i) 2 4>i°2 + 2 ( x 'l _ Xu)(x 2 - X2i)0 2 U l] + 



zi;) 2 0n (8) 

(0)l 



where "naught" superscript represent the Fermat potential at the 
position of the image and subscript represent the derivative with 
respect to the axes on the lens plane. Substitute it in equation (J7J, 
the first term after multiplication of exp[i0(xi , y;)] to the complex 
conjugate results in unity. The second term is zero from the Fermat 
principle and finally the third term as non-zero term results in mag- 
nification in the geometric optics as equation (|6j. For the third order 
4>ijk and higher orders of derivatives, the Fermat potential results 
in the optical wave feature in the light curve. 

One of the important issues in the wave optics is that, in real- 
ity, sources are not completely coherent and we may define a co- 
herent time scale of At = Aoj -1 , where Aui is the width of the 
spectrum. In this case the amplitude of the wave on the observer 
plane emitted from y is given by 



V{x,y,cf>) oc / g(uj)e 



J2R s tp(x.,y)uj 



and the magnification, taking into account the spectrum of the 
source is given by: 



2iR s tp(^,y)ui 



duj\ 



(9) 



here we set the speed of light, c = 1. For a coherent monochro- 
matic source g{ui) = ujS(u) — uio), substituting in equation Q9J, 
we recover equation Q. Assuming a non-zero temperature for a 
source, the Doppler broadening changes the spectrum of the source 
to a Gaussian distribution with the width of Aui at a given emis- 
sion line. The effect of incoherent light is diminishing wave optics 
effect, results the wave optics approaches to the geometric optics 
feature. 

An important point regarding the validity of the wave optics 
is the coherency of the light arrivi ng to the observe r plane from 
different images on the lens plane dMandzhos 198lh . In addition, 
each image should also have coherent light receiving from an ex- 
tended source. We will discuss this point in the finite size effect 
sub-section. Assuming that source is located infinitely far from the 
lens plane and has zero angular size, to examine the coherency of 
different images interfere on the source plane, we obtain the time 
delay between the two rays on the lens plane and compare it with 
the coherent time of the source. For a point source, time difference 
between the two rays received by the observer through the images 
arc 



At = 2R a [^(xn,y)-0(xi 2) y)] ; 



(10) 



where the source location is fixed and xn and x/2 are the location 
of the images. The difference between the two Fermat potentials, 
A(f> for two far images is in the order one. For the case of two 
images generated during the caustic crossing, A(f> could be in the 
order of 10~ 3 . If we quantify the time delay for the wide and close 
images in the caustic crossing, they are given as follows: 

M 



At ~ 1 x 2R S 



1(T 



M0 



(11) 
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Figure 1. The light curve of a point like coherent source with the parameters 
d = 0.8 , q = 0.1 , u = , a = , t E = 20 days , / = 1000 



and 

At ~ 0.001 x 2R S 



10" 



, M 

Me 



(12) 



We will discuss about these formulas in detail at section {3}, as 
the source gets closer to the caustic line (i.e. At — > 0). For a 
source with non-zero temperature of T, relation between the coher- 
ent time and the bandwidth of the spectrum is given by t c Auj ~ 1 
dMehta 19631) . Since the dispersion velocity of particles is related to 
the temperature by a ~ sfT and on the other hand a = Aw /u>, the 
coherent time is related to the temperature of the source as follows 
dGuenther 1990l) : 



2.8 x 10" 



fT B v 



3000 



(13) 



V(GHz) 

In the double slit experiment for a point-like source, in order to 
have diffraction pattern, the time-difference between the two light 
rays received by the observer should not be longer than the coherent 
time. Now we can constrain / = 2kR s to have diffraction features 
in the gravitational lensing, comparing equations Jilt and l l 1 2t with 
the coherent time in equation dl3t . For the wide images we have 



/ < 10' 



3000 



and for the close images we have the upper band of 



3000 



(14) 



(15) 



For the case that At > r, the results from the wave optics reduces 
to the geometric optics. A more detail analysis regarding the ef- 
fect of temporal coherency on the wave optics in the gravitational 
lensing is given in Appendix (A). 

For the case of binary system composed of a parent star and a 
planet, we can calculating equation Q. In this system, the Fermat 
potential is given by 



(x,y) 



2 (X " y) 



-m* ln(|x — x+|) — m p ln(|x — x p |), (16) 



where m* and m p are the relative mass of the star and planet nor- 
malized to the overall mass of the system and x*, x p are the loca- 
tions of lenses. We locate lenses along xi axis and put the center of 
coordinate on the center of mass of this system 
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Figure 2. The two dimensional luminosity pattern from a point source on 
the observer plane lensed by a binary system. The fringes are demonstrated 
near a critical line, while the source is crossing the caustic line. The overall 
flux would be the sum of the nearby images (as shown in this figure) plus 
flux from the incoherent wide images. 



• 0) 



,0), 



(17) 



+ ' ' * v 1 + q' 

where s± is the projected distance between the two lenses normal- 
ized to the Einstein radius and q — m p /m*. The track of source on 
the lens plane is identified with two parameters of minimum impact 
parameter, uo with respect to the center of mass and the inclina- 
tion angle of the source track a with respect to the x-axis. For the 
case of single lens, equation Q has analytic solution however for 
two point mass, we perform numerical calculation to get the light 
curve. In the next section, we also do Fourier expansion of Fermat 
potential up some relevant terms and obtain analytical wave optics 
expression. 

For a point source with a binary lens, using equation Q, we 
plot the amplification pattern as a function of time in Figure (0. 
This is a typical microlensing light curve with the wave optics fea- 
ture. Here we have two oscillating modes due to the interference be- 
tween the wide and close images. As the source gets closer to caus- 
tic line, the longer mode is magnified and after the caustic crossing 
it gets dimmer. In Figure ((2) we depict the two dimensional pattern 
on the observer plane. The relative motion of the observer with re- 
spect to this pattern produces the light curve in Figure (TJ. In the 
following section, we will consider finite coherent time and finite 
size effects of the source star in our calculations. 

2.1 Finite size effect 

The majority of microlensing observations in recent years have 
been done towards the center of Galaxy by the two main MOA 
and OGLE surveys. The main aim of these two groups is dis- 
covering the extra-solar planets by the microlensing method. The 
source for the microlensing events in the bulge due to the se- 
lection bias are relatively populated by the r ed giants compare 
to the main sequence stars dAnsari et al. 20121) . In addition to the 
visual band, red giants can also optically pump the interstellar 
medium around these stars and produce M aser emission few as- 
tronomical units far from the source star dMessineo et al. 20051 : 
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Iviemmings et al. 2005b . In contrast to the simple point like source, 
in reality source stars are extended objects. The extended sources 
not only decreases the strength of magnification in the geometric 
optics, it can also decrease enhancements of the fringes in the wave 
optics ( Schneider & Schmid-Burgk 1985b . 

The overall intensity for the exten ded source in the wave op- 
tics is given by dSchneider et al. 1992b (see also equation IA10I in 
the appendix): 



A* = 



Is My) rf y 



(18) 



where p(y) is the overall magnification of multiple images given 
by equation l |A7| > and p is the radius of source in terms of projected 
Einstein radius on the source plane. 

One of the features of the finite-size effect in the wave optics is 
that fringes with the smaller modes in the light curve dissolve as we 
increase the size of the source. This effect can be seen in the Young 
experiment when we increase the size of pinhole as the source of 
double split experiment. In the gravitational microlensing, to have 
coherent images on the lens plane, difference between the Fermat 
potential of the light rays, from the different parts of the source to 
the images on the lens plane should not be larger than light wave- 
length. Assuming L 3 as the size of source located at Di s from the 
lens plane and having an image at h off-center of the optical axis, 
the difference between the path from the two sides of the source is 



AL 



Li-L 2 = yjDf s + {h + L s ) 2 - y/Dj s + h? (19) 
hL s 
AT' 



(20) 



since the terminal point of the light rays are fixed point on the lens 
plane, the gravitational potential terms have been canceled in the 
subtraction. Difference on in the Fermat potential is related to the 
difference in the light path as AL = R s A<1>. Imposing the condi- 
tion of AL < A and substituting from equation J20l >. we obtain the 
condition of 



AO < 



Pis A 
D s h' 



(21) 



where AO is the apparent size of the source from the observer plane. 
For a single lens and high magnification events, h is in the order of 
the Einstein radius. We can write equation J2U in terms of the size 
of the source and Einstein angle as 



L s < 



Pis A 
Dot 0« 



(22) 



Using the definition of the Einstein angle, upper limit on the size of 
the source obtain as 



L s < A, 



2R^D ol 



(23) 



In order to estimate the size of source for various masses of the 
planets and wavelengths for the observation, we assuming lens is 
located at the middle of the distance between the observer and 
source where the probability of the microlensing observation is the 
highest (i.e. D a i = D; s ). For two typical planets with the masses 
of Earth and Jupiter, the Schwarzschild radius is about 1 cm and 
286 cm. Putting the source star at the Galactic bulge, D 3 — 8.5 
kpc, Table Q} shows the maximum size of the source that can pro- 
duce coherent images on the lens plane. Using A = 10 cm for the 
observation, we can observe the wave optics features for an earth 
mass planet with a typical star in the main sequence. Also we can 



A 


3 cm 


1 mm 




3.4 X 10 6 


114 


L s (Mj) 


2.0 X 10 5 


6.7 



Table 1. Maximum size of the source in kilometer produces coherent im- 
ages, located at D B = 8.5kpc far from the observer at the Galactic bulge 
observing with the two typical wavelengths of 3cm and 1mm. The calcula- 
tion is done for two planets of earth mass and Jupiter mass. 



observe a Jupiter mass planet having compact sources. In the mi- 
cron wavelengths, our limit decreases to the ten to hundred kilome- 
ter. This means that in the smaller wavelengths for observing the 
diffraction patter n, we may pro be smaller structures on source star 
such as granules dYu et al 201 lb . 

As we noted, in the binary lenses, distance between the two 
close images can be smaller than that in the single lens as the source 
approaches to a caustic line. For the wide images also , they form 
at few astronomical units far from each other. Hence at the same 
time we can have close coherent images from a planet companion 
and incoherent wide images. Hereafter, we separate magnification 
from the coherent and incoherent images into two parts of A = 



3 LIGHT CURVE NEAR CAUSTIC LINE 

Light curve of the microlensing during the caustic crossing of the 
source, can be treated by numerical and analytical methods. At the 
position of the images on the critical lines, the lens equation is valid 
hence 



4 0) - 



= 0. 



Diagonalizing Fermat potential, we can set <j>\ 2 = (P21 = 0- 
In order to satisfy singular Jacobian transformation on the criti- 
cal lines, from equation {6]l, either <pn or 4*22 should be zero. We 
set (j) 1 ^ = and tj>^ 7^ 0. Now we expand the Fermat potential 
around the critical line and ignore higher orders in x 2 as follows 
dSchneider et al. 1 992; Jaroszvski & Paczvfiski 1995b : 



(x,y) 



1 .2 
2- 



+ 



l> 222 X 2 



$1*1) + ^W-SW + <&2)Z1Z2 + ... 



Using the Fermat principle for the critical line we demand 
having singular Jacobian which implies the constrain of 

(/> l ° 2 \xi + 4> ( 2 %x 2 = 0. 

The position of the source as a function of the images can be 
obtained from the Fermat potential as follows: 



1 AO) 2 , AO) 



(25) 
(26) 



where keeping the first order term for yi and combining constrain 
from the Jacobin in j/2 term, we obtain the position of the images 
as follows: 
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(27) 



At the negative side of the caustic we have no images and at the 
positive side, two images can be formed. Now we let the source 
just cross the caustic line with an arbitrary direction. Substituting 
the position of the images in equation J24l l. we can calculate the 
Fermat potential for the nearby images during the caustic crossing. 
The difference between the Fermat potential of two images is 



a j 2 (2y 2 ) ; 

A ^3 7S 



(28) 



Here Acj> is just a function of y 2 , distance to the caustic line. How- 
ever in the observation of the caustic crossing the source can cross 
this line with an inclination, then we can replace 3/2 = sin 7 x (t — 
t c ) /tE where the time variation of A(f> is related to the transverse 
velocity and inclination angle 7 as well. Using the Fermat potential 
for a binary system in equation dl6t . we can calculate 0222 f° r me 
center of mass coordinate as follows: 



, (cm) 
w 222 



6:C2 



'p 77 



6x2 



*+xiy 



((X! - X pl ) 2 + xl) 3 



((X! - S*i) 2 + x\) 2 ((X! - X+i) 2 + x\ 



(29) 



'-'27 \ V. ) \^ 2 7 

Here x\ and X2 represent the position of critical lines where images 
form along these lines, depending on the location of the source, we 
can have double images around the parent star ( wide images) or 
around the planet (nearby images). These double images during the 
caustic crossing approach to each other. Now we need performing 
coordinate transformation of the Fermat potential from the center of 
mass coordinate to the local diagonalized coordinate at the image 
position. First we do a boost along the x-axis to one of lens po- 
sitions. The second boost is along the radial direction with size of 
Einstein radius corresponds to that lens. Finally we do rotation with 
Rij to get diagonalize </>ij . The third derivative of Fermat potential 
in the new coordinate system obtain as <p^ 2 = R2iR2jR2k<l>f : ^ > ■ 
The mentioned boosts of the coordinate to the location of the 
images around the critical line of the planet is given by x\ = 
Xpi + I Re cos 8 and X2 = R^ /Re sin 8 where R^ is the 
Einstein radius of the planet and Re is the overall Einstein radius 
and 8 is the polar angle from the position of the planet. The order of 



magnitude 0222 



at the position of the images would about 



d^rn) 

ijk 



\Rij\ ~ 1. Substituting boosts in equation l |29t , since m* ;§> m p , 
first term can dominate as R^ appears in the denominator, hence 



do) 

J 222 



,R 



E 

Re 



do) 

J 222 



(30) 



Replacing the ratio of Einstein radius of the planet to the Einstein 
radius of the star with corresponding masses and using normalized 
planet mass with m p = q/(l + q), results in 

1 

'71' 

Intuitively, comparing the double images around the planet and the 
lensing star, nearby images around the planet is more suitable for 
producing the diffraction pattern. Having smaller q makes larger 

v'222 

Now we substitute equation d30t in equation l !28t and replace 
the difference in the Fermat potential to the time difference as fol- 
lows: 



At 



2 x R s Act> 



CriTiCaUine 


y \ Caustic 
Source track 

/ B \\ 






^^^^^^^ 







5 x 10 

,t-t c 



sin 7) 



3/2/ 



t E 



-r 3/2 ( 



0.001 ' { MO> 




Figure 3. Critical lines (left panel), caustic lines (right panel) and light 
curve in the geometric optics (lower panel) for a binary system with the 
parameters of q = 0.001, Sj_ = 1. Here we take source star with the size 
of p = 0.002. Large loop around the primary lens corresponds to the par- 
ent star and smaller loop corresponds to the planet. The straight solid line 
indicates the path of the source. 



1 AOdays ' 

where (t — t c ) /1ft is the time corresponds to the relative velocity of 
the source on the lens plane normalized to one hour. We note that 
unlike to the case of single lens where At is in the order of crossing 
time of the Schwarzschild radius, in the case of bin ary lens we have 
the factor of q, decreases the corresponding time dHevl 2O10l) . On 
the other hand At approaches to zero as t — > t c . The characteristic 
time difference in the Fermat potential is about At ~ 5 x 10~ 9 
sec corresponds to l t h ~ 15 cm. This length scale corresponds 
to the frequency of 0.2 GHz. Here wavelengths larger than this 
threshold l t h, can interfere and produce the diffraction pattern from 
the images on the source plane. 

Repeating this calculation for the distant images around the 
parent star, we get a larger value for At, destructive for the wave 
optics feature. The other observable parameter which depends on 
the position of the images is the amount of the magnification during 
the lensing. Substituting equation \24\ in equation {7}, and keeping 
leading terms in the Fer mat potential, the mag nification near the 
caustic obtain as follows l lSchneider et al. 19921) : 

M~[^(f^)] 2 , (3D 
ro 

where Y = (^^)^ ~ g~ 1/6 and Ai(x) is Airy function. The 
maximum magnification also relates to the derivatives of Fermat 



! 1//3 Near the critical 
makes smaller am- 



potentialas Umax oc \<pl\'\ " I ^222 1 "'" ~ 1' 
line for the nearby images, having larger 0222 
plitude for the diffraction light curve and larger Yo produces longer 
modes. On the other hand for the distant images around the lens 
star, q is large, results in smaller 2 22 and smaller diffraction pat- 
tern with larger magnification (i.e. closer to the geometric optics). 
We remind that in our diagonalized coordinate system for the Fer- 
mat potential, the non-zero term of the second order derivative of 
the Fermat potential is 0n. On the other hand the trace of 4>ij is 
V 2 = 2 — V 2 t/>, where for two point mass lenses the second term 
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4 OBSERVATIONAL PROSPECT 
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Figure 4. The light curves correspond to the wide images (upper panel) 
and close images (lower panel) around the lens star and planet, crossing the 
caustic lines at point A and B of Figure[5] The solid line corresponds to the 
geometric optics, with the two wave optics light curves at 5 and 60 MHz 
frequencies. 



is a delta-Dirac function and can be set to zero along the caustic 
lines, far from the position of the lenses. Hence for a small distance 
from the caustic lines, <pn = 2. Therefor, the relevant parameters 
in the light curve of the wave optics are /, c/>222 and the trajectory 
of the source with respect to the caustic line. 

Figure <£3j shows the configuration of a binary lens (parent star 
with the planet) with the corresponding critical lines and caustic 
lines and the light curve in the geometric optics. Caustic crossing 
happens at two points of A and B produces nearby and distant im- 
ages. Figure 10 shows corresponding light curve in the wave optics 
and geometric optics in two wavelengths at the points A and B. 
The light curves have been obtained from the numerical Fresnel- 
Kirchhoff integration. The analytical approach in this section re- 
vealed us the physics of the wave optics in the gravitational mi- 
crolensing by a binary system. For full analysis of the light curves 
we will use the numerical calculation in the next section. 



In this section we study the possibility of follow-up observations of 
a binary microlensing system during the caustic crossing by a radio 
telescope, such as future SKA project. In a simple lens, the param- 
eters involve in the light curve are the Einstein crossing time ts, 
minimum impact parameter uo and time for the maximum magni- 
fication to. Amongst these parameters, the only parameter that has 
the physical information of the system is tE which contains the 
mass of lens, relative distance of the lens with respect to the source 
and the transverse speed of the lens with respect to our line of sight. 
Using additional information as p arallax effect due to the annual 
motion of the earth around the sun l lGould 1998) can break the de- 
generacy between the lens parameters dRahvar et al. 20031) . On the 
other hand the finite size effect of the source can also provide ex- 
tra information to break p artially the degeneracy of the problem 
dRoulet& Mollerach 19971) . 

Increasing the number of lens from one to two, increases the 
number of parameters of this system. The additional parameters 
are (i) the projected distance between the lenses normalized to the 
Einstein radius, "s±" , (ii) relative mass of the lenses, "q" and (iii) 
angle " a" defining the trajectory of the source with respect to the 
line, joining two lenses. For the binary system, we have overall six 
parameters to fit the light curve. The probability of observation of 
the binary lenses depends on the size of the caustic and for the case 
of resonance where the distance between the lenses is in the order of 
Einstein ring, we will have the maximum probability. However the 
problem with the resonance events is that despite of occasional ob- 
servation of planets, due to lack of our complete knowledge about 
the observational efficiency, it is difficult to have a statistical anal- 
ysis and extract the distribution function of the planet parameters. 
There are some efforts to have a fully-deterministic strategy by au- 
tomated se arching system to h ave the detection efficiency of the 
exoplanets dDominik et al 20100 . Having such a system to cover all 
the alerted microlensing candidates will enable us to get correct 
statistical distribution of the planet parameters. 

The other important chann el for the exoplanet obse rvation 
is the high magnification events dGriest & Safizadeh 19981) . These 
events due to the high magnification character are alerted by the 
microlensing surveys and the follow-up telescopes could monitor 
them with high sampling rate and better photometric precision. Un- 
like to the low magnification events, for these events, the detection 
efficiency function is almost known and the statistical analysis of 
the planet distribution can be applied dGould et al. 20 id) . However 
these events suffer from s o s _1 degeneracy, where we can have 
the same light curve for the close and wide caustic configurations. 
In the caustic classification of the binary systems there are three 
type of topologies for the caustics and corresponding critical lines 
dSchneider & Weifi 1986 1) , so-called the "close", "wide" and "in- 
termediate" or resonance binaries. Figure {5]l shows these three cat- 
egories of caustic lines with three different planet to the star mass 
ratio. For studying the wave optics feature during the caustic cross- 
ing, we generate synthetic light curves for the three different cate- 
gories and compare the results with the geometric optics. Our aim is 
also to study the sensitivity to wave optics signals in the resonance 
and high magnification channels. 

In order to quantify the wave optics feature in the light curve, 
we use x 2 difference from the best fit of wave optics and geometric 
optics. Assuming Oi as the error bar for each data point, [V? as the 
magnification in the geometric optics and /x' 1 " as the magnification 
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in the wave optics, the difference between the x 2 s is given by 



Ax 2 



2 

Xg.o. 



2 

Xw.o. 



E 



■/4 9J )(2/4 



(32) 



Having a threshold for A\ 2 , we can distinguish the wave op- 
tics light curve from the geometric optics. An important factor for 
determining equation ( |32t is knowing the photometric error bar 
which depends on the source flux, integration time and the size 
of the radio telescope. Amongst the various sources, red giants and 
supergiants can emit electromagnetic waves in the longer wave- 
lengths. Radio-loud Quasars at the cosmological scales are also 
bright radio sources where caustic crossing features can be ob- 
served in the cosmological scales. 

Detail study about the radio sources have been done for the 
wave optics microlensing in Heyl (2011). In what follows we 
adapt his classification. For the red giants a closer star as Arc- 
turus JPerrvman et al. 19971) shows emission at the wavelengths 
of 2cm and 6cm has 0.68 mjy and 0.28 mjy flux, respectively 
dDrake & Linskv 1986h and its spectrum is given by 



t ~ v>A I V \0.8 rkpc 2 7 

^ 24( g7^ } ( d7 } nJy - 



(33) 



Another class of stars are low mass late-type stars as asymptotic 
giants. For example, Mira is a sample of this class a nd the spectrum 
is given by jPerrvman et al. 1997llReid et al. 19971) 



U - 72 ( 



>2,kpC 2 

) ( — ) nJy. 



GHz' v D 



(34) 



Finally super giants have strong radio emission. Betelgeuse is one 
of the red super giants that is located at the distance of 197 pc 
dNewell & Hiellming 19821 ; lHarper et al. 20081) . The spectrum of 
this star normalized to the kilo-parsec distance is given by 



f" - 9 - 3 (tttft) (-75-) ^Jy 



GHz 



(35) 



The radius of super giants are in the order of a few astronomical 
unit. Assuming these stars located at the Galactic center, about 
eight kpc far from us, from equation J21l >. we can have coherent 
images for the close images on the lens plane. 

In equation d32t . we need an estimation for the photometric 
error bar. For the S KA project the nois e corresponds to Nyquist 
sampling is 0.27 Jy dSchilizzi et al. 20071 ). This sampling is defined 
according to the integration time multiplied to the band width to be 
equal to one (i.e. Av At — 1). Since noise decreases with the 
square root of time and band width as 1/ v AuAt, we can write 
the noise in terms of these two parameters as follows 

Av 



<V-aMCjgj>-"<£>-"V* 



(36) 



Now we use six parameters of the microlensing in the binary 
system as well as the source stars listed before for generating syn- 
thetic light curves with the corresponding error bars. In order to 
compare the simulated data with the geometric and wave optics 
light curves, we use the criterion of A\ 2 > 5. For an ensemble of 
light curves, we identify the caustic lines in the three categories of 
binaries as shown in Figure {3). Those light curves cross the caustic 
lines and satisfy the discrimination criterion between the geometry 
optics and wave optics are identifies in this figure. Results show 
that wave optics feature is just depends on the location of the source 
crossing the caustic lines. 

As we noted before, there are two main channels for the ex- 



< 






Figure 5. Caustic lines for three category of wide, close and intermediate 
separation of the lens from the planet with different values of q. The thick 
green area of the caustic lines produces light curves with Ax 2 > 5 from 
equation j32t . The binary lenses are indicated by two spots and scales in the 
figures are given in terms of distance between the lenses. 



oplanet observations and a qualitative study of wave optics light 
curves can be done by studying Figure (|5j. For the case of high 
magnification events, the source passes near the lens star. Having 
small q, wide and close binaries can produce almost same light 
curves and according to Figure {3J, we will not have strong wave 
optics features for this case. On the other hand for the intermediate 
regime (resonance) a larger area of the caustic lines is suitable for 
the wave optics light curve. For this case, we will have combination 
of both close and wide images on the lens plane. 

Amongst the light curves we generated in the Monte-Carlo 
simulation, a sample of light curve is shown in Figure ® with the 
cadence rate of 45 min and signal to the noise ratio of S/N — 7. 
Here the solid line indicates the geometric optics and the dotted 
line is theoretical wave optics light curves. The relevant parame- 
ters of the wave optics from equation ( 13 1 1 are the wavelength in 
the Airy function, Yq and maximum amplitude of the light, fimax ■ 
However these parameters also depend on / and 1^222 • For an en- 
semble of the light curves with uniformly distribution in the param- 
eter space, we calculate A^ 2 between the wave optics and the ge- 
ometric optics as in equation d321 >. In order study relation between 
the discriminating parameter, A\ 2 with / and ^222, we specify ar- 
eas in this parameter space, satisfying A^ 2 > 10 with the source 
sizes of p — 0.001, 0.005 and 0.002, in Figure {JJ. Upper part of 
these curves are excluded. Having smaller / means longer wave- 
length for the observation of small mass lens. On the other hand 
larger 1^222 is also relates to the wave optics from the time differ- 
ence between the images. Smaller size sources satisfying the spatial 
coherency are suitable for the wave optics microlensing. Knowing 
the transit time scale of the fringes changes with respect to the ob- 
server At, we can determine Yo(f, 2 2 2 ) = Ar/t_E. On the other 
hand yL max {4>222) can be measured directly from the light curve. 

The physical parameters involve in the wave optics light 
curves are the the overall mass of the system Alt, q and the angle 
of the source trajectory, a. On the other hand form the observations 
in the visual telescopes, we can find the best light curve fit with six 
parameters of the theoretical model and possible degrees of degen- 
eracy. Having extra information from the wave optics will constrain 
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Figure 6. A part of theoretical light curve from the geometric and wave 
optics. The data points are simulated according to the amount of noise from 
SKA. Here we take -£ = 7. The parameters of the lens is similar to that 
in Figure (4) with the observational frequency of 10 GHz frequency and 45 
minute integration time. 



100000 



1000 















p=0.002 




p=0.005 




p=0.01 



20 40 60 80 100 120 

4*222 

Figure 7. Area above the curves in parameter space of (/ — ^222) space 
has been excluded for having wave optics features with the criterion of 
Ax 2 > 10. The size of the source star normalized to the Einstein radius is 
mentioned in the legend of Figure. Smaller stars are favorable fort the wave 
optics microlensing. 



Mt, q space and subsequently we can have better confidence from 
measurement of s±. 

Now we want to examine the sensitivity of the wave optics 
signals to the main physical parameters of the binary system, q and 
s±. In our Monte-Carlo simulation, we can identify percentage of 
the events that satisfy the condition of A\ 2 > 10. According to 
Figure 1J8}, the fraction of events is given in terms of percentage 
and suitable area of the parameter space for the wave optics signals 
is in the resonance region where the separation between the lens 
and planet is around the Einstein radius. The physical reason as 
we discussed in Figure {5), is that in the intermediate regime, the 
close images favorable for the wave optics features can be formed. 
In order to estimate the overall number of the microlensing events 




Projected distance 

Figure 8. Efficiency function in s, q space that can produce wave optics fea- 
tures on the observer plane. The criterion for the wave optics from equation 
(32} is taken A\ 2 > 10. The simulation is done for the source star with the 
finite size of p = 0.005 and observation at the frequency of v = 10GHz. 
The six parameters of the lens are generated with the uniform distribution. 
Red lines separate the parameter space according to the definition of the 
wide, close and intermediate binaries from Erdl & Schneider (1993). 



with the wave optics signal, we should multiply this histogram to 
the diagram contains the abundance of the binary events. 

Finally we want to extract the physical information from a 
typical wave optics light curve, assuming that observation is done 
both in the visual and radio wavelengths. From the visual wave 
length observation, fitting with the geometric optics we can extract 
q and s± and trajectory of the source. On the other hand from the 
wave optics our relevant variable are / and 1^222 ~ 1/ -Jq. 

In what follows we assume a synthetic data taken by the ra- 
dio observation along with the geometric observation in the visual 
band. Here parameters of the lens as the apparent distance between 
the lenses as well as the mass ratio and trajectory of the source can 
be obtain from the geometric optics, however there is a degree of 
degeneracy in the best fit to the data. From the wave optics the de- 
gree of freedom appear in the light curve is / which relates to the 
total mass of the lenses. This parameter can be fixed by the wave 
optics observation and subsequently the total mass of the system 
can be identified. In Figure l[9} we fit the synthetic data with the 
theoretical wave optics light curve. Here the theoretical value of / 
is 49475 and assuming observation at v — 10 GHz, from the like- 
lihood function we get / = 51400^ ]^qq. Using k in the definition 
of /, we can extract the overall mass of the binary system. 



5 CONCLUSION 

In this work we study the wave optics effect of the gravitational mi- 
crolensing by a binary system composed of lens star and compan- 
ion planet. Having planet in the lensing system causes generating 
close images during the caustic crossing, a suitable configuration 
of the images on the lens plane to produce diffraction pattern. This 
effect is in analogy to the double slit Young experiment in the as- 
tronomical scales. We derived the light curve for the binary system 
in the framework of the wave optics showing that it depends on the 
derivatives of the Fermat potential and /. 

Taking red giants and supergiants as the source stars of the 
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Figure 9. Synthetic light is generated by the cadence of 30 min and photo- 
metric precision of 0.15, fitted with the wave optics theoretical light curve. 
The initial value of / is taken 49475, using v = 10GHz observational 
frequency. The likelihood function shows the best fir with la and 2a level 
of confidences. With la, we obtain / = 51400^ 



-,+1635 
-1600 



gravitational microlensing in the galactic scales where there is a 
natural preference for observing these type of sources in the mi- 
crolensing, we suggestion using SKA future project for the ob- 
servation of wave optics signals in the light curve. In this obser- 
vational program, radio observations can accompanying the mi- 
crolensing follow-up telescopes in the visual bands towards the 
center of Galaxy. These two observations in the long and the shorter 
wave lengths can be a complimentary observations of the binary 
microlensing events, result is breaking the possible degeneracy of 
binary system. Amongst two observational channels for the exo- 
planet observation by the microlensing so called (i) high magni- 
fication channel and (ii) resonance channel, we showed that wave 
optics observations is applicable for the resonance events where a 
larger area of the caustic lines is favorable for manifestation of the 
wave optics microlensing. Finally we discuss on sample of syn- 
thetic microlensing light curves and application of wave optics and 
geometric optics to identify the parameters of the lens with the bet- 
ter accuracy. 
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APPENDIX A: COHERENCY IN THE WAVE OPTICS: 
GRAVITATIONAL LENSING 
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In order to show the effect of coherent time on the wave optics 
in the astronomical scale, we compare the time delay between the 
light rays arriving to the observer with the source coherent time. We 
adapt the wave optics formalism in Schneider et al. (1992). Starting 
from the amplitude of the electromagnetic field on the lens plane, 
it can be obtained on the observer plane as 



o if<j>(x,y) , 2 



(Al) 



Doing the Taylor expansion of the Fermat potential around the im- 
age at the stationary point and diagonalizing the Fermat potential 
results in 



The integration of equation l lA2t is 
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where det | J\ is the determinant of <j>ij in equation (6), in the di- and simplifying this equation results in 
agonalized coordinate and n refers to the type of images and can „ 

be n = 0, 1, 2, depending on the number of focal points transverse fi — — J fi(s)ds (A10) 
from the source to the observer ( I Arnold 19891) . Now if we have N 
images from the lensing, the overall amplitude would be 

V = ^±Y, (A4) 

We note that this equation is valid while the source is out of the 
caustic lines (i.e. det \Ji\ ^ 0). If we take into account the spec- 
trum for the source, replacing f with 2ujR s , the overall amplitude 
is given by 

N . 

The overall magnification is given by equation Q, hence the over- 
all magnification after averaging over the time is given by 

N 

fj, = y^^i (A6) 

/ f I ^2 i[2fl s ^(0 (O) -0 (O) )-f (n,- n ,)] , 

The argument in the phase is f((f>i — </4°') = ^At, assuming 
a Gaussian spectrum for the source, for the case that the temporal 
coherent time is shorter than the time difference between the light 
rays, the cross terms will be zero and we reach to the geometric op- 
tics. For the simple case of two images, N = 2, putting a Gaussian 
spectrum for g(ui), we obtain the following result for the magnifi- 
cation 

(i = /tii + ii 2 + 2e y T> y/juj^ cos[cj At - — (m - n 2 )],(A7) 

here we replace the width of the spectrum with the coherent time 
as follows Aoj = -, For At <C r, we have an oscillating mode, 
reveals the wave optics feature from the superposition of the waves. 
In contrast, for At ^> r the exponential term vanishes the oscillat- 
ing term and we will have geometric pattern for the magnification. 

In the case of extended source, each point on the source has 
contribution on every point of the observer plane. Hence we for this 
case we can write equation l |A5t for each point of source assigning 
it by V(s). The overall amplitude can be written as 



= / / V , dSdS ' x (A8) 



i| K s ( B ^( 8 )- U '^( S '))-K-n 3 >/2] < > ^du*. 

where S is the area of the source and averaging is taken over the 
time, hence < g(cj)g(u)') >= g(cj) 2 S(uj — a/). Since the differ- 
ential elements on the source are spatially uncorrelated, the cross 
terms in equation l lA8t from the different elements will be canceled 
and only light rays propagating from one element of the source con- 
tribute in this summation. Mathematically we may consider 

<4°>( a )^V) >~ 6(s-s'). 
Hence equation jAS\ can be simplified to 

i r r N i[it»w(0f ) W-^° ) (s))-(n i -n 3 )7r/2] 

** = c / ds I XI , , T == 9{u) 2 du, (A9) 
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